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In this short survey paper, we shall discuss certain recent results in classical gravity. 
Our main attention will be restricted to two topics in which we have been involved: the 
positive mass conjecture and its extensions to the case with horizons, including the Penrose 
conjecture (Part I), and the interaction of gravity with other force fields and quantum- 
mechanical particles (Part II). 

I. Positive Mass Conjecture and Related Topics 

One of the most difficult problems in classical relativity is to understand how and when 
singularities form. In the 1960's, Hawking and Penrose proved that the existence of a 
closed trapped surface in an asymptotically flat spacelike hypersurface gives rise to a 
singularity in space-time. However, no proof based on pure partial differential equation 
arguments was found, and many questions remain unanswered. 

Given an initial data set (gij,Pij) on a three-dimensional manifold so that gij is asymp- 
totically Euclidean and pij (the induced second fundamental form in an embedding) falls 
off asymptotically, it is interesting to ask the following questions. 

1. When will such an initial data set contain a closed trapped surface? If so, how to 
locate it? 

2. If the initial data contains no closed trapped surface, how to tell whether such a 
surface will appear at a later time under the evolution of Einstein's equations? 

3. If we assume that the trace of p^ is zero and that under the evolution this condition 
is preserved, will a singularity occur without the existence of a closed trapped surface 
for all time? 

4. If a singularity does occur, what is the structure of the null geodesies in a neigh- 
borhood of the singularity, and what is the structure of the curvature tensor in this 
neighborhood? What is the criterium on the initial data set for the curvature to 
blow up at the singularity? 

5. Can one define physically relevant local (or quasi-local) quantities such as mass 
and angular momentum to describe regions in a strongly gravitationally interacting 
space-time? For example, when two bodies interact, what is the binding energy 
and what is the mass of the resulting configuration? How can one estimate the 
gravitational radiation for strongly interacting bodies? How can one justify the 
linearized theory of gravitational radiation? 
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For all the above questions related to singularity formation, one usually studies only generic 
initial data. However, it has been a difficult problem in nonlinear partial differential 
equations to understand how to perturb away the singularity. 

For all these questions, it would be good if the known class of spherically symmetric 
solutions of the Einstein equations were rich. Except for the Schwarzschild case, such 
solutions cannot be vacuum solutions. Hence to consider these questions, one is forced to 
couple gravity to other matter fields. For the case of a massless scalar field, Christodoulou 
|Il] , |l~3ll ) has studied the question of the formation of singularities quite exten- 
sively. If a singularity exists, it is located at the origin. While much is known in this case, 
the details of how the singularity forms is still poorly understood. (When naked singu- 
larities form, one would like to know the behavior of null geodesies.) Based on numerical 
studies, Choptuik Q has found the new phenomena that, for a one-parameter family of 
initial data, the mass function exhibits some critical phenomena similar to those which 
occur in statistical mechanics, at the time when a black hole forms. However, a detailed 
theoretical study is lacking. An interesting consequence of the above study is that after 
gravitational radiation, the space-time is either time-asymptotic to the flat space-time 
or to the Schwarzschild space-time. This raises an interesting question when we couple 
gravity to a Yang-Mills field or to Dirac spinors: what is the possible asymptotic state of 
spherically symmetric initial data? Would those stable coupled solutions found by us (cf. 
Section II below) be the only possible states? When we want to extend the spherically 
symmetric case to an axisymmetric geometry, the space-time is far more complicated. 
While it is clear that angular momentum may be used to make many configurations sta- 
ble, the number of degrees of freedom is large and it is difficult to find solutions of gravity 
coupled to other fields. (For stationary black holes with a vacuum background, it has to 
be the Kerr solution.) It is still not known whether one can find multi-black holes which 
can be stabilized by the addition of angular momentum. 

Beyond axisymmetric solutions, Bartnik [Q, ||] proposed a class of initial data sets which 
can be foliated by round spheres. Using this ansatz, he was able to parametrize a large 
set of initial data having zero or nonnegative scalar curvature. (For initial data set, if it is 
a maximal slice, the scalar curvature is always nonnegative.) According to his numerical 
study, this ansatz has been very useful in understanding radiation of a single black hole. 
Perhaps the theoretical study for critical data in this class would be interesting. 

Let us now turn to general space-time with no spherical symmetry. We restrict our- 
selves to asymptotically flat space-times. In this case, we have asymptotic space-time 
Lorentzian symmetry. Based on this asymptotic symmetry, it is well known that one can 
define the concept of mass and linear momentum associated to each initial data set (which 
is invariant under the Lorentzian symmetry at asymptotic infinity ||). About twenty years 
ago, Schoen-Yau |33| (subsequently |36| and others) proved the positive mass conjecture 
which says that the total (mass, linear momentum) is a non-spacelike four-vector. The 
total mass is therefore always nonnegative. It is zero only when the space-time is flat. 

The positivity of the mass says that the trivial space-time is stable (the dynamic 
stability among a class of reasonable initial data has recently been demonstrated by 
Christodoulou and Klainerman [Q|). However, the nonlinear stability of the Schwarzschild 
solution is still unknown. Based on the "Cosmic Censorship conjecture," Penrose proposed 
an inequality relating the total mass of the black hole to the area of the outermost horizon. 
It says that among all initial data sets with fixed mass, the time-symmetric Schwarzschild 
solution initial data set has the largest area for its outermost apparent horizon. 

In general, if the initial data set is a maximal slice for the space-time, the scalar 
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curvature of the three-dimensional manifold is non-negative. In such a case, the conjecture 
of Penrose was recently settled by Huisken-Ilmanen [25], obtaining the optimal result only 



under the assumption that the outermost black hole is connected. It was based on an idea 
of Geroch that the Hawking (quasi-local) mass is monotonic along an evolution of a surfaces 
T t which starts from the hole to the sphere at infinity. The evolution is governed by the 
requirement that it moves the surfaces along the normal direction and with magnitude 
minus the inverse of mean curvature. Geroch notices that for the sphere at inifinity, the 
Hawking mass is simply the total mass of the initial data set while at the black hole, the 
Hawking mass is, up to a universal constant, the square root of the area of the black hole. 
Hence if the flow of the surface Tt exists, the Penrose conjecture would then be proved. 
Huisken and Ilmanen developed the mathematical framework in which these ideas could 
be made precise. However, the flow exhibits jump phenomena and much care is needed 
to assure that the inequality jumps in the right manner. Much more recently, H. Bray Q 
has been able to improve the result in the case of a non-connected outermost horizon by 
a new method, partly relying on the ideas of Schoen-Yau [33] and the curvature estimates 

§• 

For the proof of the Penrose conjecture, one still must answer the question as to when 
the initial data set is not maximal. It would also be nice to see the corresponding inequality 
for the Bondi mass (total mass after radiation). 

Besides the total mass and linear momentum, an important conserved quantity is 
angular momentum. This was studied extensively by Ashtekar ||. One needs to study 
the relation between angular momentum and other conserved quantities such as the total 
mass and linear momentum. It seems reasonable to believe that the total mass should 
dominate the square of the angular momentum if the initial data set is nonsingular. 

To better understand angular momentum, Huisken and Yau |2(| defined the concept 
of center of mass of an initial data set. It is Lorentzian invariant and, remarkably, under 
the evolution of Einstein's equations, the velocity of the center of gravity is the linear 
momentum divided by 2m, where m is the total mass of the initial data set. 

One hopes to study all possible naturally conserved quantities and the relations among 
these conserved quantities, when the initial data set is non-singular. It is always inter- 
esting to know how radiation effects all those quantities. For an isolated gravitational 
system, what would be its asymptotic state after radiation? We conjecture that the time- 
asymptotic state is just the superposition of several known stationary solutions including 
the charged Kerr black holes and the static coupled solutions found above (e.g. when we 
are coupling the Einstein equations to the Yang-Mills, Dirac particles, or a real scalar 
field). 

The global behavior of the Einstein system is difficult to study, partially because 
we do not have (quasi-)local quantities which behave well under time evolution. The 
Hawking mass is one such example. It is monotonic in some directions. Unfortunately, it 
is not positive in general. For certain important closed surfaces, which are obtained by 
minimizing area under a volume constraint, Christodoulou and Yau were able to prove 
the positivity of the Hawking mass [15]. However, they assumed that the scalar curvature 



is non-negative. It would be nice to remove this assumption. 

If one considers sufficient conditions for the formation of black holes in a general set- 
ting, the best theorem is due to Schoen-Yau |Q. This says that, by suitably defining the 
diameter d(£l) of a region $7, then if the matter density in the region Q is greater than 
d~ 2 up to a universal constant, a closed trapped surface can be found. This implies that 
a black-hole type singularity exists. In this theorem, the existence of black holes results 
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from the condensation only of matter. It would be desirable to include the contribution 
of gravitation effects. Namely, it is interesting that in the argument by Shoen-Yau, only 
the lower bound of the first eigenvalue of the operator — A + | R is used. In the time 
symmetric case (i.e. with py =0), \ R is the local matter density. It would be nice to 
see if this method can be extended to the general case, in the sense that the spectrum of 
some operator can be used to yield a condition for the formation of black holes. 

II. The Interaction of Gravity with Other Force Fields and 
Dirac Particles 



According to Einstein's Theory of General Relativity, gravity is described geometrically 
through Einstein's equations. The understanding of gravity has been driven by the dis- 
covery of special solutions of these equations. The most important examples are the 
Schwarzschild solution, the Kerr-Newman solution, and the Friedmann-Robertson- Walker 
solution [|J. Particularly interesting effects are obtained when one couples gravity, as ex- 
pressed through Einstein's equations, to other fundamental force fields. The simplest such 
example is the Reissner-Nordstrom solution resulting from the coupling of gravity to elec- 
tromagnetism (Maxwell's equations). This solution, like the Schwarzschild solution, has 
an essential singularity at the origin. The generalization to non-abelian Yang-Mills fields 
led to the discovery of Bartnik and McKinnon (BM) § (see also [27, ^0[) of everywhere 
regular solutions. This came as a surprise because several results for related systems led to 
the conjecture that such solutions cannot exist. Indeed, neither the vacuum Einstein equa- 
tions, nor the pure Yang/Mills equations have non-trivial static, globally defined, regular 



solutions [17]. The existence of these solutions depends on the coupling of the different 



force fields, whereby the attractive gravitational force is balanced by the YM repulsive 
force. But this balance is rather delicate; for example, the BM solutions are unstable with 
respect to small perturbations [j35|. Other interesting solutions of Einstein's equations 
result from coupling gravity to quantum mechanical matter fields. The case of a complex 
scalar field was considered by T.D. Lee et al [28], who found soliton-like solutions modeling 
(bosonic) stars; see too D. Christodoulou [11] who studied the gravitational collapse of a 
real massless scalar field. 

We report here on recent work (see []T8[|— ]2T|] for details) of a different type of coupling; 
namely, gravity coupled to both quantum mechanical particles with spin (Dirac parti- 
cles), and to an electromagnetic field. We first study the resulting Einstein-Dirac-Maxwell 
(EDM) equations for a static, spherically symmetric system of two fermions in a singlet 
spinor state. We find stable soliton-like solutions, and we discuss their properties for dif- 
ferent values of the electromagnetic coupling constant. We note too that the inclusion of 
gravity has a regularizing effect on solutions, in the sense that our solutions are more reg- 
ular than one would expect from a naive analysis of the Feynman diagrams; see [^]. We 
then study black- hole solutions for these equations (see |27], |3l[]), and we find, surprisingly, 
that under rather weak regularity conditions on the form of the event horizon, the only 
black-hole solutions of the EDM equations are the Reissner-Nordstrom (RN) solutions. 
That is, the spinors must vanish identically. Applying this to the gravitational collapse to 
a black hole of a "cloud" of relativistic spin-^-particles, our result indicates that the Dirac 
particles must eventually disappear inside the event horizon. We also show that the Dirac 
equation has no normalizable, time-periodic solutions in a RN black-hole background. The 
physical interpretation of this result is that the Dirac particles cannot remain on a periodic 
orbit around the black-hole. This result has recently been extended to an axisymmetric 
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black hole geometry [p^ j. 

In our study of the coupled EDM equations, we employ a special ansatz for the spinors. 
In this ansatz, we do not assume that the Dirac particles are in a spherically symmetric 
state; indeed they are allowed to have angular momentum. However, we arrange (2j + 1) 
of these particles in such a way that the total system is static and spherically symmetric. 
(In the language of atomic physics, we consider the completely filled shell of states with 
angular momentum j. Classically, this multi-particle system can be thought of as several 
Dirac particles rotating around a common center such that this angular momentum adds 
up to zero.) Since the system of fermions is spherically symmetric, we obtain a consistent 
set of equations if we also assume spherical symmetry for the gravitational and electric 
fields. We can thus separate out the angular dependence, and the problem then reduces 
to a system of nonlinear ODEs. 



1 The EDM Equations 

The general Einstein-Dirac-Maxwell equations are 
R) - \R5) = -8itT; , (G-m)tf a = , V k F^ k = 4vre^^G% a , (1) 

where Tj is the sum of the energy-momentum tensor of the Dirac particles and the Maxwell 
stress-energy tensor. The G J are the Dirac matrices which are related to the Lorentzian 
metric via the anti-commutation relations 

gi k (x)l = ^{G j (x),G k (x)} = ^(G j G k + G k G j )(x) . 

F^ k denotes the electromagnetic field tensor, and ^> a are the wave functions of fermions 
of mass m and charge e. The Dirac operator is denoted by G, and it depends on both the 
gravitational and electromagnetic field; for details see [18, |19| ]. 

We now specialize to the case of static, spherically symmetric solutions of the EDM 
system ([[]). In polar coordinates (i, r, i?, (p), we write the metric in the form 

ds 2 = -L= - J- dr* - r 2 {d$ 2 + sin 2 dtp 2 ) (2) 
1 [r) z Ayr) 

with positive functions T and A. Depending on whether we consider particlelike solutions 
or black-hole solutions, the region of space-time which we consider is r > 0, or r > f > 0, 
respectively; in the latter case, we assume that r = f is the event horizon. We always 
consider solutions for which the metric (|2|) is Minkowskian, 

lim A(r) = 1 = lim T(r) , (3) 

r— >oo ) — >oo 

and has finite (ADM) mass; i.e., 

lim - (1 — Air)) = p < oo . (4) 

r — >oo 2 

In the static case, the fermions only generate an electric field, and thus we may assume 
that the electromagnetic potential A has the form A = (—0,0), where (ft = <ft{r) is the 
Coulomb potential. 
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The Dirac operator G can be written as 
= iG*M + B(x) 



where 7*, 7 r , 7 , and 7^ are the 7-matrices in polar coordinates, in Minkowski space; 
namely 



7 

7* 
1* 



7 

7 1 cos "& + 7 2 sint? cos <p + 7 3 sini? s'mip 
- ^— 7 1 sini? + 7 2 cost? cos ip + 7 s cos 1 !? sin ip 

(— 7 2 sin 93 + 7 3 cos ip 

rsinu V 



where 






-a 1 



a' 




1,2,3, 



and a 1 denote the Pauli matrices. 

In analogy with the central force problem in Minkowski space [29], this Dirac operator 
commutes with: a) the time translation operator idt, b) the total angular momentum 
operator J 2 , c) the z component of the total angular momentum J z , and d) with the 
operator 7°P, where P is the parity. Since these operators also commute with each other, 
any solution of the Dirac equation can be written as a linear combination of solutions 
which are simultaneous eigenstates of these operators. We use this "eigenvector basis" to 
separate out both the angular and time dependence, and to calculate the total current 
and energy momentum tensor of the Dirac particles. Using the ansatz in [jl8| , [19 , [20[| , we 
can describe the Dirac spinors using two real functions a, /3. We arrive at the following 
system of ordinary differential equations for the 5 real functions a, /?, A, T, and <j>: 



A a' 

~A0 

r A' 
V 

2rA — 
T 



r 2 A< 



,. 2j + 1 
± a 



((w - e<f>)T + m)P 
2j + l 



2r 



p 



2r 

((u! — e<p)T — m) a =F 

AT 

l-A- 2(2j + l)(u - e(p)T 2 (a 2 + p 2 
A - 1 - 2(2 j + l)(w - ecp)T 2 (a 2 + f3 2 
+2{2j + 1) mT (a 2 



± 2 



r 2 AT 2 \<P'\ 2 
(2j + l) 2 



Tap 



p 2 ) + r 2 AT 2 



>\2 



-(2j + l)e{a 2 + p 2 



V 



2rA + r 2 A 1 A' 



(6) 

(7) 
(8) 

(9) 
(10) 



Equations (§) and (fjj) are the Dirac equations (the ± signs correspond to the two possible 
eigenvalues of 7°-P); (||) and (0) are the Einstein equations, while Maxwell's equations 
reduce to the single equation ([Uj) . Here j = \ , | , . . . , the constant w enters via the plane 
wave dependence of the spinors; namely exp(— iut), and as for the general equations (||), 
m and e denote the mass and charge, respectively, of the fermions. We also require that, 
in addition to (Q), (0j), the electromagnetic potential vanishes at infinity, 



lim ( 

r— >oo 







(11) 
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Since the equations (|^)~([l0|) are invariant under the gauge transformations 

4>(r) — > (f>(r) + k , u> — > iv + en , kGIR, (12) 

we see that (0) can be fulfilled by a suitable gauge transformation, provided that (j) has 
a limit at infinity. 

In Sections [l|Q, we shall be concerned with two different types of solutions of equations 
(|(f)-(|i~0|); namely particlelike solutions (smooth solutions defined for all r > 0), and black 
hole solutions (solutions defined for all r > f > 0, where A{f) = and A(r) > for all 
r > f; r = r is the event horizon). In the first case, we require the following normalization 
condition on the spinors: 

/ (a 2 + B 2 ) —= dr = 1 , (particlelike) (13) 

Jo V A 

while in the second case we require that for all tq > f, 

< / (a 2 + f3 2 ) —= dr < oo (black holes). (14) 

Jr V A 

These conditions are necessary in order that the Dirac spinors define physically meaningful 
wave functions. 



2 Particlelike Solutions 

In this section we shall describe our numerical construction of particlelike solutions for 
the equations (|6|)-(|l0|). For simplicity we shall restrict ourselves to the case j = 1/2. We 
shall also discuss the stability and properties of the ground state solutions for different 
values of the electromagnetic coupling constant (e/m) 2 . We shall show that solutions exist 
even when the em coupling is so strong that the total interaction is repulsive in the non- 
relativistic limit. In addition, for small em coupling, (e/m) 2 < 1, we shall show that stable 
particlelike solutions exist for small values of m, and using certain topological techniques, 
we show that this stable solution becomes unstable as m increases. 

The construction of particlelike solutions is obtained via a rescaling argument (see 
g|, OH)- The idea is to weaken the conditions (|), (0), and ([0| to 

^ (a 2 + B 2 ) —= dr < oo , ^ lim T(r) < oo , lim Mr) < oo , 
Jo \J A r ^°° r ^°° 

(15) 

and instead set 

T(0) = 1 , 0(0) = , m = 1 . (16) 

This enables us to use a Taylor expansion around r = 0, and we obtain the following 
expansions near r = 0: 

a(r) = air + 0{r 2 ) , 3(r) = 0{r 2 ) 

A{r) = l + 0(r 2 ) , T{r) = 1 + 0{r 2 ) , <f>(r) = 0{r 2 ) . 

Solutions to our equations now depend on the three real parameters e, u>, and a%. For 
a given value of these parameters, we can construct initial data at r = 0, and using the 
standard Mathematica ODE solver, we "shoot" for numerical solutions of the modified 
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system (|6|)-([lO|), (|i~6|). By varying uj (for fixed e and a%), we can arrange that the spinors 
(a, j3) tend to the origin for large r, and the conditions (|j) and (|l5|) also hold. 

Given a solution (a, (3, A, T, <p) of this modified system, we consider the scaled functions 

a(r) = ^a(Xr) , 0{r) = ^(Ar) 

A(r) = A(Xr) , f(r) = r~ l T{\r) , 0(r) = r 0(Ar) . 

By direct computation, these functions satisfy the original equations and the 

equations (|3|), (||), and (|l3|), provided that the physical parameters are transformed ac- 
cording to 

m = Am , w = Arw , e = Ae , (17) 

where the scale factors A and r are given by 

A = (47t/ (a 2 + I3 2 ) —= dr) , r = lim T(r) . 

V Jo vA J r ^°° 

Finally, the condition (jll|) can be fulfilled by a suitable gauge transformation. Notice that 
the parameter (e/ra) 2 = e 2 is invariant under the above scaling. It is thus convenient 
to table (e/m) 2 (and not e 2 ) as the parameter used to describe the strength of the em 
coupling. We point out that the above scaling technique is used only to simplify the 
numerics; for the physical interpretation, however, we must always work with the scaled 
(tilde) solutions. Since the transformation from the un-tilde to the tilde variables is one- 
to-one, our scaling method yields all the solutions of the original system. From now on, 
we shall only consider the scaled solutions, and for simplicity in notation, we shall omit 
the tilde. 



3 Properties of the Particlelike Solutions 

We have found solutions having different rotation numbers n = 0,1,2,... of the vec- 
tor (a, 0). In the non-relativistic limit, n is the number of zeros of the corresponding 
Schrodinger wave functions, and thus n = corresponds to the ground state, n = 1 to 
the first excited state, and so on. However due to the nonlinearity of our equations, n no 
longer has this simple interpretation. For simplicity in what follows, we shall only discuss 
the n = solutions. The graphs of a typical such solution is shown in Figure |]. For each 
solution, the spinors (a, (3) decay exponentially to zero at infinity. We interpret this to 
mean that the fermions have a high probability to be confined to a neighborhood of the 
origin. In view of this rapid decay of the spinors, our solutions asymptotically go over into 
the spherically symmetric RN solutions of the Einstein-Maxwell equations Jl]], as r — > oo. 
That is, for large r, 

Mr) ~T"<r)~l-*Z + W . 

In other words, our solution, for large r, looks like the gravitational and electrostatic field 
generated by a point particle at the origin with mass m and charge 2e. Note that in 
contrast to the RN solution, however, our solutions have no event horizon or singularities. 
One can understand this from the fact that we consider here quantum mechanical particles, 
rather than point particles. Therefore the wave functions are de-localized according to the 
Heisenberg Uncertainty Principle, and so the distributions of matter and charge are also 
de-localized, thereby preventing the metric from forming singularities. In general, we can 
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parametrize solutions by the rest mass m, and the energy uj of the fermions. In Figure |2[ 
we plot the binding energy m — oj versus m for different values of the parameter (e/m) 2 , 
and we see that m — u is always positive, indicating that the fermions are in a bound 
state. For weak em coupling, (e/m) 2 < 1, the curve is a spiral which starts at the origin. 
The binding energy decreases for fixed m and increasing (e/m) 2 , since the em repulsion 
weakens the binding. The mass energy spectrum when (e/m) 2 <C 1 becomes similar to 
the case of the Einstein-Dirac equations (without the em interaction); see J18| |. We can 
use linearization techniques to show numerically that for small m, if (e/m) 2 < 1, the 
solutions are stable with respect to spherically symmetric perturbations. For larger values 
of m, we can investigate the stability using Conley index theory (see |J2|), where m is 
taken to be the bifurcation parameter. This technique shows that the stability/instability 
of a solution remains unchanged if m is varied continuously and no bifurcations occur. 
Moreover, at bifurcation points, the Conley index theory provides a powerful technique to 
analyze changes of stability. Using this, we find that all solutions on the "lower branch" of 
the spiral curves A and B of Figure ^ (i.e., on the curve from the origin up to the maximal 
value of m), are stable, and all solutions on the "upper branch" are unstable. 

From Figure ^, we see that this form of the mass energy spectrum changes when 
(e/m) 2 ~ 1, the regime where, in the classical limit, the electrostatic and gravitational 
forces balance each other. To better understand this situation, we take the non-relativistic 
limit in our EDM equations. To do this, we fix (e/m) 2 , and assume that e and m are 
small. In this limit, the coupling of the spinors to both the gravitational and em forces 
becomes weak: A, T ~ 1 and ~ 0. The Dirac equations imply that to ~ m and a ^> (3. 
Thus the EDM equations go over to the Schrodinger equation with the Newtonian and 
Coulomb potentials; namely, 

( — — A + e<f> + mV) ^ = E ^ (18) 
\ 2m / 

-AU = -8vrm|^| 2 , -A(f> = 8vre|^| 2 , (19) 

where E = uj — m, 'l'(r) = a(r)/r, V(r) = 1 — T(r), and A is the radial Laplacian on 
H 3 . From tfHfy we see that the Newtonian and Coulomb potentials are multiples of each 
other; namely V = — m/e (j). Thus if (e/m) 2 > 1, the total interaction is repulsive so that 
the Schrodinger equation ( |l8| ) has no bound states. It follows that in the limit of small 
m, the EDM equations have no particlelike solutions, if (e/m) 2 > 1. This means that the 
mass-energy curves in Figure || can only start at m = if (e/m) < 1. This is confirmed by 
the numerics (Figure curves C, D, and E). For (e/m) 2 = 1, the curve tends tom-w = 
as m — > oo. 

If (e/m) 2 > 1, Figure § shows that the EDM equations admit solutions only if m is 
sufficiently large, and smaller than some threshold value where the binding energy of the 
fermions goes to zero. 

We can also consider the total binding energy p — 2m, where p is defined in (Q). In 
Figure || we plot p — 2m versus m, for various values of (e/m) 2 . If (e/m) 2 < 1, p — 2m 
is negative for the stable solutions, while p — 2m > if (e/m) 2 > 1. This indicates that 
if (e/m) 2 > 1, such solutions should be unstable because energy is gained by breaking up 
the binding. 
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4 Non-Existence of Black Hole Solutions 



As we have noted in the last section, particlelike solutions of the EDM equations in a given 
state (e.g. the ground state) cease to exist if the rest mass m of the fermions exceeds a 
certain threshold value m s . The most natural physical interpretation of this statement is 
that if m > m s , the gravitational interaction becomes so strong that a black hole would 
form. This suggests that there should be black hole solutions of the EDM equations for 
large fermion masses. In this section, we shall show that this intuitive picture of black 
hole formation is incorrect. In fact, we prove that under weak regularity conditions on the 
form of the horizon, any black hole solution of the EDM equations must either be the RN 
solution (in which case the Dirac wave function is identically zero) , or the event horizon has 
the same general form as the extreme RN metric. In the latter case, we show numerically 
that the Dirac wave functions cannot be normalized. It follows that the EDM system does 
not admit black hole solutions. Thus the study of black holes in the presence of Dirac 
spinors leads to unexpected physical effects. If we apply this result to the gravitational 
collapse of a "cloud" of Dirac particles, our result indicates that the Dirac particles must 
eventually disappear inside the event horizon. 

In order to establish these results, we first recall what is meant by black hole solutions 
of the EDM equations. These are solutions of Equations (|6|)-(|To|) defined in the region 
r > f > 0, which are asymptotically flat (so that (|3|) holds), and have finite (ADM) mass 
(so that @) holds), and satisfy the normalization condition (|l4|). In addition, we assume 
that A(r) > for r > f, and lim r \ f A(r) = 0, while T{r) > and lim r \f T(r) = oo. 

We make the following three assumptions on the regularity of the functions A, T, and 
(j) on the form of the event horizon r = r: 

(I) The volume element ^/|det gij\ = r 2 A~2 T _1 is smooth and non-zero on the horizon; 

i.e., 

T- 2 A-\ T 2 A G C°°([f, oo)) 

(II) The electromagnetic field tensor is = d{Aj — djAi\ we assume that the strength 

of the em field tensor F i5 F ij = -2|^| 2 A T 2 is bounded near the horizon. In view 
of (I), this means that we assume 

\<t>'(r)\ < c\ , f<r<f + e 

for some positive constants c%, £ > 0. 

(III) The function A(r) obeys a power law; i.e. 

A(r) = c(r-r) s + 0({r - r) s+1 ) , r>r (20) 

for some positive constants c and s. 

A brief discussion of these assumptions is in order. Thus, if (I) or (II) were violated, then an 
observer freely falling into a black hole would feel strong forces when crossing the horizon. 
Assumption (III) is a technical condition which seems sufficiently general to include all 
physically relevant horizons; for example s = 1 corresponds to the Schwarzschild horizon, 
and s = 2 corresponds to the extreme RN horizon. However, (III) does not seem to be 
essential for our non-existence results, and with more mathematical effort, we believe that 
it could be weakened or even omitted completely. 
Here is the main result in this section. 
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Theorem 4.1 Any black hole solution of the EDM equations ffi)-(lC) which satisfies the 
regularity conditions (I) -(III) either is a non-extreme RN solution with a(r) = = j3{r), 
or s = 2 and the following expansions are valid near the event horizon r = f: 



A(r 
T(r 

<fi(r 

a(r 
0{r 



A (r - rf + 0((r - f) 3 ) 
To(r-f)- 1 + 0((r-f)°) 

- + <Po (r-f) + C((r-f) 2 ) 
e 

a {r-ff + 0{(r-f) K+1 ) 
Po(r-f) K + 0((r-r) K+1 ) , 



(21) 
(22) 

(23) 

(24) 
(25) 



with positive constants An, Tq, and real parameters (fin, an, and (3q. The exponent k 
satisfies the constraint 



1 

2 < K 



A -lJ m 2_ e 2 (/) 2 T 2 + 



(2j + l 

V 2f 



and the spinor coefficients an and (3q are related by 

2j + l 



a 



1 An K ± 



2f 



/?o (m - e^ r ) , 



(26) 



(27) 



where rb ' refers to the two choices of the signs in f|qV-f|l6|j. 

We shall now outline a proof of this result; we first consider the case that the exponent 
s < 2 in ©. 

Lemma 4.2 Assume that s < 2 and i/iai (a, f3,A,T,(fi) is a black-hole solution where 
(a, /?) ^ 0. T/ien i/iere are constants c, e > satisfying 

1 



c < a(r) / + y S(r)^ < 
Proof: According to (§) and @, we have 



^A±{a 2 + f3 2 ) = 2 
ar 



r < r < r + e 



(28) 




2j + l 



2r 
-m 



< 



,2\ 5 




(c^ + /^) 



(29) 



The uniqueness theorem for ODEs implies that (a 2 + /3 2 )(r) > for all r, r < r < r + e, 
for any e > 0. Dividing ([29] ) by \/3 (a 2 + /3 2 ) and integrating from r > f to f + e gives 



log((a 2 + /3 2 )(f + £ )) - log((a 2 + /? 2 )(r)) 



A - 2(t) Urn 2 + 



(2j + 1 
t 2 



\2\ 2 



(30) 



Since s < 2, ( p0| ) implies that i 2 is integrable on r < r < r + e, so that the integral in 
(]30"|) is majorized by 



r+£ 



A 4m 2 + 



(2j + 1) 
t 2 



1 

2\ 2 
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and this yields (|28| 



We can now dispose of the case < s < 2; namely, we have 

Proposition 4.3 If < s < 2, then the only black hole solutions of the system 
are the non-extreme Reissner- Nordstrom solutions. 

Proof: We assume that we have a solution such that (a,(3)(r) ^ 0, and show that this 
gives a contradiction. 

The last lemma implies that the spinors are bounded near r = f. From @ and (||), 
we find 

d (2? + Fl 2 

r-f(ylT 2 ) = -4(2j + 1) (u - e<P) T 4 (a 2 + (3 2 ) ±2 ' T 3 a(3 

dr r 

+2(2j + l)mT 3 (a 2 - P 2 ) . (31) 

Assumption (I) implies that the left side of ([3TJ) is regular so the same is true of the right 
side. Since T — * oo as r V f, we conclude that 



lim (w-ed>(r)) = . (32) 



From Maxwell's equation 



^" = 'J (« 2 + /? 2 ) " V 2 ^T]' ^ , (33) 

we see that (I) implies that the coefficient of (j)' is smooth. If s > 1, vl _1 is not integrable 
at f, so that \(j)'\ is unbounded at f, thereby contradicting (II). Thus s < 1, and integrating 
( |33"1) twice and using (^) gives near r = f the following expansions: 

0'(r) = ci (r - r)~ s+1 + c 2 + C((r - f)~ s+2 ) , and 

0(r) = Cl (r - f)~ s+2 + c 2 (r - f) + - + £)((r - f)- s+3 ) . 

e 

Using these in (||), and noting that A and r 2 AT 2 \(j)'\ 2 are bounded near r = f, and that 
(oo — e4>) = 0(r — f), and T 2 (a 2 + (3 2 ) ~ (r — f) _s , s < 1, we see that the rhs of (§) 
is bounded near r = f. On the other hand, the lhs of (|8|) diverges near r = f since 
ryl'(r) = (r — f) _s+1 ; this contradiction completes the proof. ■ 



In the case s > 2, we first prove the following two facts (cf. 

lim(r-f)~3 (a 2 + (3 2 ) = (34) 

r\r 

and 



i 



lim |<£' (r) | = f _1 lim^~2T _1 > . (35) 

From (|3"5|), we find that 

(w — ecj))(r) = c + d (r — f) + o(r — f) , 

where d = e/f lim r \fA~2 T _1 > 0. Thus (u — e<j))T diverges monotonically. From (^) 
and (0), this implies that liminf r \ f (a 2 + f3 2 ) > 0, thereby contradicting (34). Thus if 
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s > 2, there are no solutions of (|^)-(|To|). 

Proof of Theorem £A: We must only consider the case that s = 2 and (|2l|), ( p2|) hold. 
From (34) we see that lim r \f a 2 + (3 2 = 0, and we can show that {u — e<f))T cannot 
diverge monotonically near r = r (see ^lj). But ( j35[) shows that (a; — e<f>) has a Taylor 
expansion near r = f with a non-zero linear term. Thus (^) holds, the constant term in 
the Taylor expansion of (u> — e<p) vanishes, and lim r \ f (a; — e4>)T = A, where from (|3F 
|A| = f _1 lim r 
variable 



A 2 T 1 > 0. As in [21 1, we may write the Dirac equations in the 



u r 



r ln(r 



and apply the stable manifold theorem to conclude that a and (5 satisfy the power laws 
and (H) yields that K > \. Using (|2l])-(|25]) into (§ and (0) gives 



' Aq k ao 



± 1 «o + (e^olb - m) /? 
zr 

2j + 1 

-(e^oTo + m)a =F n _ /?o 
2r 



which are equivalent to (|2~6|) and (p7|). This completes the proof of Theorem 4.1 



Notice that in the case of non-zero spinors (s = 2), Theorem 4.1 places severe con- 
straints on the behavior of black hole solutions near the event horizon, in the sense that 
since k > ^, the spinors decay so fast at r = f, that both the metric and the em field behave 
like the extreme RN solution on the event horizon. Physically speaking, this restriction 
to the extremal case means that the electric charge of the black hole is so large that the 
electric repulsion balances the gravitational attraction, and prevents the Dirac particles 
from "falling into" the black hole. Of course, this is not the physical situation that one 
expects in the gravitational collapse of, say, a star. However, extreme RN black holes are 
physically important since they have zero temperature [24], and can be considered to be 
the asymptotic states of black holes emitting Hawking radiation. It is thus interesting to 
see if the expansions (21)-(^5|) yield global black hole solutions of the EDM equations. 

This question is especially interesting since in the next section we shall show that 
for an extreme RN background field, spinors satisfying the expansions (^4|), ( p5| ) cannot 
be normalized. The question thus becomes whether the influence of the spinors on the 
gravitational and em field can yield black hole solutions with normalized spinors. This is 
a very difficult question because one must analyze the global behavior of these solutions 
of the EDM equations. Our numerical investigations show that the answer to the above 
question is negative; namely solutions either develop a singularity for some r > f, or the 
spinors (a, j3) are not normalizable. We thus conclude that the expansions (|2T|)-(f25|) do 
not give normalizable solutions of the EDM equations. 



5 Dirac Particles in a Reissner-Nordstrom Background 

In this section, we shall consider solutions of the EDM equations where we fix the back- 
ground metric and em field to be an RN solution. Near a collapsing black hole one might 
guess that Dirac particles might get into a static or time periodic state. However, we 
shall show that, in contrast to the classical situation, the Dirac equations do not admit 
any normalizable time-periodic solutions; in particular, they admit no normalizable static 
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solutions. We do not assume any spatial symmetry on the wave functions. This result can 
be physically interpreted as saying that Dirac particles can either disappear into the black 
hole of escape to infinity, but they cannot remain on a periodic orbit around the black 
hole. We note that it is essential for our arguments that the particles have spin. In fact, 
in the case where the particles do not have spin, the Dirac equation must be replaced by 
the Klein- Gordon equation, and our arguments fail; c.f. pgfl. 
The RN metric can be written in polar coordinates as 

ds 2 = (l - y + ^ dt 2 - (l - *L + ^ dr 2 - r 2 (d$ 2 + sin 2 dip 2 ) , (36) 

where p is the (ADM) mass of the black hole, and q its charge. The em potential is of the 
form (—0,0) with Coulomb potential 

0(r) = q - . (37) 



In the "non-extremal" case (q < p), the metric coefficient (1 — £ + %) vanishes twice, 
and thus there are two horizons < vq < r. If q = p, the metric is called an extreme 
Reissner-Nordstrom (ERN) metric and has a single horizon at r = p. If q > p, the above 
metric coefficient is non-vanishing, and so the metric does not describe a black hole; this 
case will not be considered. 

We consider time-periodic solutions, noting that static solutions are a special case. 
Since the phase of the Dirac wave function \& has no physical significance, we define \& to 
be periodic with period T if for some real f2, 

*(t + T, r, 0, tp) = e~ inT #(i,r, *!?,</?) . (38) 

Our main theorem in this section is the following: 

Theorem 5.1 i) In a non-extreme RN background, there are no normalizable, time- 
periodic solutions of the Dirac equation, ii) In an ERN background, every normalizable, 
time-periodic solution of the Dirac equation is identically zero in the region r > p. 

We shall begin by deriving conditions which relate the wave function ^ on both sides 
of the event horizon. We first consider the case of a non-extreme RN background, and 
analyze the behavior of near the event horizon. For this, we begin by studying the 
behavior of ^ in a Schwarzschild background metric, and we shall also consider the Dirac 
equation in different coordinate systems. This is done with the aim of passing to Kruskal 
coordinates, in order to remove the "Schwarzschild singularity." 

The Schwarzschild metric is 



ds 2 



l-^E)dt 2 - (l-^) 'dr 2 - r 2 (dtf 2 + sin 2 tfV 



r / 



where p is the (ADM) mass, and the event horizon is at r = 2p. Some straightforward 
calculations (see pl]| ) shows that outside the horizon (r > p), the Dirac operator can be 
written as 

Gout = | y | + y ( iS » + i ( S - 1) + t s) + rf± + .7^ , (39) 
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where 

S(r) = 

r 

The normalization integral is considered over the hypersurface t = const; i.e. 

(^|*)* ut := ( {^m.^S-U^x , (40) 
Jm 3 \b 2p 

where B2 P denotes the ball of radius 2p about the origin, and ^ = ^*7° is the adjoint 
spinor. In the region r < 2p, the Dirac operator is given by 

with corresponding normalization integral 

/ (Vj r y)(t,x) S- 1 d 3 x . (41) 

JB 2 p 

Our description of spinors in this coordinate system poses certain difficulties. Namely, 
since the t variable is space-like inside the horizon, the normalization integral (^l|) is not 
definite since the integrand is not positive. Thus we can no longer interpret the integrand 
as a probability density. Moreover, the Dirac equations corresponding to the operators 
Gi n and G ou t describe the wave functions inside and outside the horizon, respectively. But 
it is not evident how to match the wave functions on the horizon. To handle these issues, 
we remove the singularity at r = 2p by going over to Kruskal coordinates. Recall (see |l]D 
that Kruskal coordinates u and v are defined by 

^ 1 e^p cosh [ — j for r > 2p 

» = < ' *V ^ )\ P > (42) 
1 e 4 p sinh I — I for r < 2p 



2p V4p 



T JL. ( t 

1 e 4 P ginh — ) f or r > 2p 

2p UpJ 

r~ j_ ft 

1 e 4 p cosh — ] for r < 2p 

2p \4p 1 



(43) 



The horizon r = 2p maps to the origin u = = v, and the singularity r = maps to the 
hyperbola v 2 — u 2 = 1, v > 0. In Kruskal coordinates, the metric fl36| ) becomes 

ds 2 = r 2 (dv 2 -du 2 ) - r 2 (d-d 2 + sin 2 dip 2 ) 

where f~ 2 = e 2 ? . Taking v and u as time and space variables, respectively, and 
noting that the metric is regular at the origin, we can extend the Dirac operator smoothly 
across the origin. A straightforward computation gives the Dirac operator in Kruskal 
coordinates as 

G = l t (fi-i- + -f (d v r) - l - d v f] + y f/i A + I (/ ( du r) - 1) - I d u f 



dv r 2 ) \ du 

+ 17% + 17^ • (44) 
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Observe that the Dirac operator is smooth across the event horizon. Moreover, the nor- 
malization integrals (40) and ( fH] ) on the surface t = become 



Jn 



m 
where 

n = {u = 0, < v < 1} U {v = 0, u > 0} , 
is the normal to H pointing into the region u > 0, v > 0, and G 3 are the Dirac matrices 

G v = / 7 * , G u = ff , G & = 7^ , G v = 7^ . 

We remark that for smooth solutions of the Dirac equation, one can use current conser- 
vation 

Vjf(?$ = , (45) 

to continuously deform the hypersurface H keeping fixed the value of the normalization 
integral. For example, one can deform TL to Ti as depicted in Figure |3], thereby avoid- 
ing integrating across the horizon. On the other hand, one must exercise extreme care 
whenever a solution of the Dirac equation is singular near the origin. 

As shown in [21], the Dirac operator in Kruskal coordinates can be written as 

G = U G out U- 1 = U G in U- 1 , (46) 

where U is the time-dependent matrix 

U(t) = cosh (J^j 1 + sinh (J-^j 7 * Y , (47) 

and the Dirac operators G ou t and G m in Kruskal coordinates are 

-S+ l -S')^ - V + <7 * * +*y* * . 
r 2 J r oi) dip 

It follows that the Dirac operators G ou t and G; n can be identified with the Dirac operator 
G in the region 

K = {u + v > 0, v 2 -u 2 < 1} . 

We next see how solutions of the Dirac equation inside and outside the horizon match 
on the horizon, u = = v. To do this, we first study the behavior of these solutions on 
the horizon. Let us first consider static solutions of the Dirac equation, so 

tf(t s r,0,p) = e-*"*tf(r,0,y) • 

We assume that ^ is a solution of the Dirac equations (Gj n — m) = and (G ou t — m) = 0, 
and that \P is smooth on both sides of the horizon r < 2p and r > 2p. Using (f46|) and 
((47|), we have 

V(u,v,0,(p) = U(t)e~ iujt y{r,$,v) , 
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where r and t are determined implicitly from u and v in the usual way (see Q). This 
implies that ^ is only defined in 1Z, and solves there the Dirac equation (G — m)\P = 0. 
Since we are only considering black holes, we demand that VP vanishes in the half-plane 
u + v < 0; thus we must analyze solutions of the form 



U(t) e~ iujt *(r,tf 5 ¥>) for u + v > 0, u / t> 
for u + i> < . 



Such a wave function might be singular along the lines u = ±v, in which case VP must 
satisfy the Dirac equation in a generalized sense. An analysis carried out in [^] shows 
that \& must satisfy the two matching conditions 



lim(7* + 7 r ) \e\i^(t,2p + e,$,<p) 







\e\i{*(t,2p + e,0,<p) - V(t,2p-e,0,tp)) 
= o(l + |e|4 ^(t,2p + e,tf,tp)) as e -> 0. 



(48) 



(49) 



Note that since these only depend on the local behavior of Vl/ near the horizon, they are 
also applicable when we are in the case of a non-extreme RN background having event 
horizon at r = 2p. 

We now consider Dirac particles in a RN background. Since the gravitational and EM 
background fields are spherically symmetric and time independent, we can separate out 
the angular and time dependence of the wave functions via spherical harmonics and plane 
waves in the usual manner and, as shown in 21], we obtain the following two component 
Dirac equations: In regions where the t-variable is time-like, 



dr i ku > 



(u — ecj>) 



S 



2j + l 
2r 




jku) 



(50) 



and in the regions where t is space-like, 



S dr® jkuJ 



[CO 




+ i 



m 



jkuj 



(51) 



In these equations, 



S(r) 



1 



^ 4. t 



j = |, |, . . ., k = —j, —j + 1, . . . ,j, and the ± signs correspond as before to the two 
eigenvalues of the operator 7°P (cf. Section |l]). Here we have chosen for the Dirac wave 
functions the two ansatz' 



(52) 



jkui 



-iuit 



-iuit 



X* + x ^i(r) 



(53) 
(54) 
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with 2-spinors ^ jkuJ , and X +1 are denned by 



k j + l-k fc -l / 1 \ j + l + fc fc+l / \ 

" + '+i 1 o J " V~*7+2~ '+i ^ i J • (56) 

where YJ fc are the usual spherical harmonics, / = 0, 1, 2, . . ., /c = — I, 

We shall show that the matching conditions (f48|),(|49D do not yield normalizable, time- 
periodic solutions of the Dirac equation. This will be done by showing that, for every 
non-zero solution of Dirac's equation, the normalization integral outside and away from 
the horizons 

| *)* = / "^7*^ S- 1 d 3 x , (57) 

JM 3 \B2 ri 

is infinite for some t. Note that for a normalizable wave function, this integral is the 
probability that the particle lies outside the ball of radius ri, and thus cannot exceed 1. 
So if (57) is infinite, the wave function cannot be normalized. 

Now assume that ^ is a T-periodic solution of Dirac's equation. Expanding the periodic 
function e~ lQt ^(t, r, if) in a Fourier series, and using the basis (||),(||) yields 

*(i,r,tf )¥ >) = Yl , (58) 

n,j,k,s 

where s = ±, and u(n) = f2 + 2irn/T. Using the orthonormality of the spinors \ • i 1 ; 

the 

J ± 2 



integral (|5q) becomes 

\ 2r i n,n' j,k,s 

In order to eliminate the oscillating time dependence of the integrand, we average over 
one period (0, T) to get 

IjVl*)^ = J2 (*5m»)I*5* W («)) • 

" n,j,k,s 

For a normalizable wave function, this expression is finite, and hence all summands are 
finite: For all s, j, k, n, 

(*ifc W (n) I *ifc W (n)) < 00 • ( 59 ) 

We shall show that (|5^) cannot hold for non-trivial solutions of the Dirac equation; for 
this we begin with 

Lemma 5.2 The function l^jj^COl 2 has finite boundary values on both horizons, and if 
it is zero on one horizon, then it is identically zero. 

Proof: For simplicity, we omit the indices j, k, and u. Choose 5, < 6 < ro, and notice 
that the t-direction is time-like on the intervals (5, ro) and (ri,oo). In these regions, we 
can use (|^) to obtain 

si$±|2 (r ) = <s A$±, $±> + <$±, s 4U±> 

dr dr dr 



(|$±| 2 -|$±| 2 ) - 4m Re ((<!>*)*$*) , 
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so that 

_ c |$±|2 < S^-l^l 2 < C\^\ 2 

dr 

where c = 2m + (2j + 1)/S. Dividing by \^^\ 2 and integrating gives, for 5 < r < r' < 7*0, 
or 7*1 < r < r', 

-c f S" 1 < logl^l 2 T < c f r S- 1 . (60) 

J r r J r 

In the region tq < r < ri, (51) gives similarly 

S ^.|$±|2( r ) = <S^-*±, cD±> + «D± tfiU±> = , 
dr dr dr 



since the square bracket in ( |5l|) is an anti-Hermitian matrix. Thus l^ 1 *^ 2 is constant in this 
region, so (|60[ ) trivially holds for 7*0 < r < r' < r\. Since S^ 1 is integrable on the event 
horizons, (|60| ) shows that \^^\ 2 has finite boundary values on each side of the horizon, and 
these are non-zero unless if vanishes identically on the corresponding region (<5, ro), 
(r ,ri), or (7*1,00). 

We now use ( |53|) and ([54] ) in the matching condition (|49| ) to get for j = 0, 1, 



$ ± (rj +e) -$ ± (r i -e) = o(l + |$ ± (r i + e)|) , e -> 0. 

Since ^^(r)! 2 has 2-sided limits as r^, we conclude that these limits must coincide at rj; 



i.e. 



lim |$ ± (r,-+e)| 2 = lim |$ ± (r 1 - - e)| 2 . 

Using (|6^) again, we conclude that the wave function vanishes on the entire interval (5, 00), 
if it is zero on rj . This completes the proof since 5 was arbitrary. ■ 



The final step is to use current conservation (cf. pi 

V i ¥G i * = (61) 

to study the decay of ^jku n \( r ) at infinity, and to prove Part i) of Theorem 5.1 : 

Theorem 5.3 (radial flux argument) Either vanishes identically, or the normal- 
ization condition Km) is violated. 



Proof: For simplicity, we again omit the indices s, j, k, and u. Suppose that VP ^ 0. 
For 7*1 < r < R and T > 0, let V be the annulus outside the horizon r, given by 
V = (0, T) x {B 2R \ Bar). Using ph, we find 



= / V,(^G%) W^j d 4 x 

= / cftr 2 S(r) / (^ 7 r ^)(t,r) - / dt R 2 S(R) f f* 7 r #) (*, B) 
Jo Js 2 Jo Js 2 

p2R p __ 

- / dss 2 5 _1 (s) / (*7**)(t,r) 

J2r JS 2 



t=T 
t=0 



Since the integrand is static, the last term vanishes, and we conclude that the radial flux 
is independent of the radius, 

r 2 Sir) ( (tfyWr) = R 2 S(R) [ f^ 7 r *)(B) . (62) 
Js 2 Js 2 



19 



Using (|53|) and (pjf ), we have 

r 2 5(r) / (^7 r ^)(r) 



g 2 



S 2 



**(r) 



i 
-i 



The matching condition (5q), expressed in terms of <3? gives 



$(r) 



(63) 



lim 

n<r— »n \ i —1 



Using this, we have from 



lim r 2 5(r) / (f 7 r ^)(r) 



ri<r-^ri 



S 2 



lim 



lim 



1 t 

-i 1 

1 
-1 



1 i 
i -1 



l$l 



lim / |$| 2 / 

p i<i — ►n J52 



since $ is finite and non-zero on the horizon T\. 

Now we consider the radial flux for large R. Since the flux is independent of R, we 
have from the last inequality 

< lim \R 2 S(R) [ C$j r m(R)\ 

lim \R 2 S~ l {R) I (*7**)(i2)| , 



R-*oo 



S 2 



because the metric is asymptotically Minkowskian. Thus the integrand of our normaliza- 
tion integral 



O I *) 



2ri 



dRR 2 S- 1 (R) / (*7*^)(i?) 



s 2 



converges to a positive number, so that the normalization integral is infinite. 



We have thus proved Part i) of Theorem 5.1. For Part ii), the case of an extreme RN 
background field, we use a quite different method; c.f. |2l|| . 

We remark that, using Chandrasekhar's separation method, the results in this section 
can be extended to the axisymmetric case. Namely, for a quite general class of axisymmet- 
ric black-hole geometries, including the non-extreme Kerr-Newman solution, it is proved 
in [23] that the Dirac equation admits no normalizable, time-periodic solutions. 
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Figure 1: Solution of the EDM equations for parameter values (e/m) 2 = 0.7162, m = 
0.7639, u = 0.6807, p = 1.15416 (a'(0) = 0.05361). 
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Figure 3: Total Binding Energy p- 2m for (e/m) 2 = (A), 0.7162 (B), 0.9748 (C), 1 (D), 
and 1.0313 (E). 




Figure 4: Kruskal Coordinates 
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